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ABSTRACT 

Constraints on cosmological parameters from upcoming measurements with the Mileura Widefield 
Array—Low Frequency Demonstrator (MWA-LFD) of the redshifted 21 cm power spectrum are fore¬ 
cast assuming a flat ACDM cosmology and that the reionization of neutral hydrogen in the inter- 
galactic medium occurs below a redshift of z = 8. We find that observations with the MWA-LFD 
cannot constrain the underlying cosmology in this scenario. In principle, a similar experiment with 
a 10-fold increase in collecting area could provide useful constraints on the slope of the inflationary 
power spectrum, n s , and the running of the spectral index, a Sl but these constraints are subject to 
the caveat that even a small reionization contribution could be confused with the cosmological signal. 

In addition to the redshifted 21 cm signal, we include two nuisance components in our analysis related 
to the systematics and astrophysical foregrounds present in low-frequency radio observations. These 
components are found to be well separated from the signal and contribute little uncertainty (<30%) 
to the measured values of the cosmological model parameters. 

Subject headings: Cosmology: Early Universe, Galaxies: Intergalactic Medium, Radio Lines: General, 
Techniques: Interferometric 


1. INTRODUCTION 

Hydrogen gas in the intergalactic medium (IGM) is a 
promising tool for studying some of the most interesting 
topics in astrophysics, including the nature of the first 
luminous objects, the processes of structure formation, 
and cosmology. Prior to the epoch of reionization (EOR), 
when radiation from the first luminous sources reionized 
the IGM, the 21 cm hyperhne line from neutral hydrogen 
is predicted to be visible against the cosmic microwave 
background (CMB) with relative brightness temperature 
of order ±10 mK. At high redshifts, the spin temperature 
of HI is expected to be less than the CMB temperature, 
making the line visible in absorption, whereas at lower 
redshifts, after radiation from the first luminous objects 
has begun to heat the IGM but before it has reionized 
it, the line should be evident as emission. 


As the primordial hydrogen cools following recom¬ 
bination and later reheats, density contrasts are re¬ 
vealed as fluctuations in the bri ghtness tempera ture of 


the redshifted 21 cm line (Sunvaev & Zclclovich 

1972 

Hoga 

n & Reed 1979tlScott & Rees 1990!: Iliev et al. 

2002 

2003; 

Loeb & Zaldarriaga 2004|; Barkana & Loeb 2005b) 


At high redshifts, the power spectrum of these fluctua¬ 
tions is expected to follow closely the dark matter power 
spectrum and measurements from this epoch could, in¬ 
dependently from CMB and galaxy clustering experi¬ 
ments, constrain cosmological models. At lower red¬ 
shifts, as the first luminous objects ionize their surround¬ 
ings, voids are expected to appear in the diffuse emission 
(iMadau et al.lfl997t iTozzi et al.ll200d: ICiardi &; Madaul 

I200.1t IZaldarriaga et all12004 iFurlanetto et al.l l2004f l. 
Measurements of the power spectrum during this epoch 
would help reveal the characteristics of the first luminous 
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objects. 

Recent efforts have distinguished theoretical aspects 
of redshifted 21 cm observations that are particu- 
larly relevant for con straining cosmological models. 
iBarkana fc Loebl ( 2005cf ) discuss the ability of redshifted 
21 cm observations to constrain the geometry of the high 
redshi ft universe be twee n recombination and rei oniza¬ 
tion: lAli et al.l (|2005l ) and IBarkana fc Loebl (l2005ah have 
shown that differences in the line-of-sight versus angu¬ 
lar components of a redshifted 21 cm power spectrum 
measurement can be used to separate primordial density 
perturbations from features caused by t he radiati ve pro- 
cesses responsible for reionization; and IBarkana i 200(» : 
has discussed the application of the Alcock-Paczynski 
(AP ) test to redshifted 21 cm measurements. Addition¬ 
ally, fBarkana^Toefcj (j2005bi ) consider the effects of the 
earliest galaxies on the r edshifted 21 cm fluctuations and 
iNaoz &; Barkanal (l2005l l discuss using redshifted 21 cm 
observations to study the thermal history of hydrogen 
gas by detecting a cutoff in the power spectrum due to 
thermal broadening of the line. 

The experimental challenges of detecting and char¬ 
acterizing neutral hydrogen in the IGM at high red¬ 
shift are significant. While reaching the sensitivity limit 
needed to detect the redshifted 21 cm radio background is 
readily achievable by very small experiments, statistical 
measurements of the fluctuation power spectrum require 
the collecting area of the much larger arrays—MWA- 
LFD, LOFAR, and PA ST—currently under development 
(jBowman et al.ll200fih . Even after the required collecting 
areas have been reached, anticipated challenges remain. 
Astrophysical foreground contaminants are five orders of 
magnitude brighter than the redshifted 21 cm emission. 
Removing the foreground signatures from observations 
requires extremely precise observational calibration. 

In this paper, we expand on earlier calculations of the 
sensitivity of the first generation arrays to forecast their 
ability to place constraints on cosmological models from 
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measurements of the redshifted 21 cm power spectrum. 
In order to treat the best-case scenario and to simplify 
the analysis, we have chosen to assume that reionization 
has not yet begun by the target redshift range of 8 < 
2 < 10. This assumption represents the most optimistic 
scenario generally consistent with the existing evidence 
about when reioniz ation began from quasar absorption 
line measurements dDior govski et al .11200 ll: iBecker et al.1 
120011 : iFan et all 120031 iWvithe fc Loebl 120041) and from 
Tho mson scatteri ng measurements by the WMAP satel¬ 
lite (ISpereel et all 120061) . These measurements suggest 
reionization occurred at redshifts z > 6 and z ~ 10 , re¬ 
spectively. 

We begin in Section 2 by reviewing the observational 
process of measuring the highly redshifted 21 cm power 
spectrum and describing the fiducial experiments. In 
Section 3, the method to forecast constraints on cos¬ 
mological parameters is described in terms of statistical 
errors using a Fisher matrix treatment of the full three- 
dimensional power spectrum. The calculation marginal¬ 
izes over several cosmological parameters and two an¬ 
ticipated contributions related to the astrophysical fore¬ 
ground contaminants. The results are discussed in Sec¬ 
tion 4 and encapsulated in Figure [2j 

A similar study ha s been performed in parallel by 
iMcQuinn et al.l 1 2006( 1. While we analyze the character¬ 
istics relevant to constraining cosmological models with 
the initial generation of experiments and in the presence 
of two nuisance contributions, their analysis explores the 
benefit of combining redshifted 21 cm measurements with 
information from other experiments, such as WMAP and 
Planck. Together, the two efforts provide a thorough 
overview of the potential of future redshifted 21 cm power 
spectrum measurements to contribute to cosmology. 


2. THE MEASUREMENT 

The basis of the statistical measurement of the red¬ 
shifted 21 cm power spectrum using low-frequency, wide- 
held radio observations has been d ev eloped in the lit¬ 
eratur e bv iMorales k, Hewittl (I2004T). IZaldarriaea et al.l 
(l2004f h iMoralesI ~i 20051) . and iBowman et al.l (|2006l . 
These efforts are built on the similar approach 
employed for interfer ometric me asurements o f CMB 
aniso t ropies (IWhite et all 119991 : iHobson fe Maisingerl 
12002 ;: iMvers et al.l l2003f T~Here. we review the relevant 
properties of the measurement; we direct the reader to 
the previous works for additional details. 

Neutral hydrogen is optically thin to the 21 cm line. 
Thus, through the combination of their angular and 
spectral responses, the visibility measurements of low- 
frequency radio arrays inherently sample the emission 
from a three-dimensional volu me of space at high red- 
shift (|Morales &; Hewitt] 120041) . The measurements may 
be represented in a fully Fourier domain and expressed 
in coordinates convenient for studying cosmology (k = 
ki,k 2 ,ks) or in coordinates more directly related to in¬ 
strumental considerations ( u,v,p ). 

The fluctuation power spectrum of redshifted 21 cm 
emission in the sampled volume of space is mapped to 
an instrumental response by the convolution of the power 
spectrum, Phi{ k), with the instrumental window func¬ 


tion, W(k), according to (|Morales fe Hewittl I2004D 

C HI ( k) = <|A/(k)| 2 > = J P HI (k)\W(k-k')\ 2 d 3 k\ 

(!) 

where the window function is given by the Fourier trans¬ 
form of the instruments angular and frequency response, 
and is very sharply peaked for the the cases we will con¬ 
sider (see Section 2.2). 

There is an inherent uncertainty in the measurement of 
the redshifted 21 cm power spectrum from cosmic sample 
variance. The uncertainty can be estimated by dividing 
the observed three-dimensional Fourier space into a large 
number of independent cells, where the volume of each 
cell is approximately the size of the window function, 
Wik). For the regime of interest in our calculations, 
the variance in the measured power spectrum is then 
CYj{ k) ~ C HI (k)5ij , where the indices i and j run over 
all the independent cells in the sampled volume. 


2 . 1 . Contributions from Astrophysical Foregrounds 


The redshifted 21 cm emission is not the only source 
of power in low-frequency radio observations. Astro- 
physical foregrounds are several orders of magnitude 
stronger and dominate the measured signal. The fore¬ 
grounds include free-free and synchrotron emission from 
the Galaxy, extragalactic point sources, and free-free 
emission from elections in the IGM. 

The bright foregrounds manifest themselves in several 
ways. First, the Galactic emission dominates the ther¬ 
mal noise in the measurements, especially at lower fre¬ 
quencies. Although much of this large-scale emission is 
expected to be resolved-out by interferometric observa¬ 
tions, its effect on antenna temperature remains. 

The contribution to the measured power spectrum due 
to thermal noise is substantial, but white ((C^k)) = 
constant ), and therefore should be readily removed. In 
our calculations, we will assume that it has been removed 
(but allow for an imperfect subtraction by including an 
offset parameter in the model). After this subtraction, 
the thermal uncertainty per independent cell in the three- 
dimensional measured power spectrum can be approxi- 
mated, in instrumental coordinates, using (lMoralesll2005l . 
His Eqn. 11) 


[ C Ui ( U )] 


= 2 


/ 2kBT sys 
\ edAdp 


1 


B n(u) t 


( 2 ) 


where dA is the physical antenna area, dp is the inverse 
of the total bandwidth, ks is Boltzmann’s constant, T sys 
is the total system temperature, e is the efficiency, B is 
the total bandwidth, n is the time average number of 
baselines in an observing cell, and t is the total obser¬ 
vation time. This equation is derived from the corre¬ 
sponding relationship for interfer omet ric m easurements 
of CMB anisotropies in IWhite et al. (1999. Their Eqn. 
16). Although the measured power spectrum in this case 
is three-dimensional, the parameters in Equation [2] are 
approximately independent of frequency within the ob¬ 
serving band. The thermal uncertainty per independent 
cell, therefore, is taken also to be independent of p (and 
£ 3 ). Additionally, for the MWA-LFD, the time-averaged 
visibility distributions are expected to be dense and to 
have nearly circular symmetry, thus n(u) depends only 
on \Ju 2 + v 2 . 
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The second important effect of astrophysical fore¬ 
grounds is that they produce their own signatures in the 
measured power spectrum. The dominate foreground 
signature is expected to be due to extragalactic (and 
Galactic) point sources, which introduces a distinct com¬ 
ponent in the observations that is very bright (rang¬ 
ing up to ~1000 K for individual sources) and highly 
structured on the sky. Fortunately, these foreground 
contaminants are expected to be distinct from the red- 
shifted 21 cm signal in their angular correlation functions 
and in their frequency behavior; and recent studies indi¬ 
cate that statistical techniques should provi de meth ods 
to se p arate them from the d e sired signal (lOh k, Mackl 
20031 iDi Matteo et al.l 2004 IZaldarriaea et al.l 2004 : 


Morales fe Hewitt! 12004 Santos et al.l 120051 : IWang et alJ 
20061 iMorales et al J l2005l~ 


We will assume that the foregrounds have been cleaned 
from the measured signal. Even the best subtraction al¬ 
gorithms leave residual traces of the contaminate, how¬ 
ever, and a perfect model may produce an imperfect sub¬ 
traction simply due to thermal uncertainty in the mea¬ 
surement. Such a subtraction error would result in a 
power spectrum signature followi ng the shapes of the 
model components (jMorales et al.l i2005h . The principle 
residual signature due to this kind of subtraction error 
for extragalactic point sources can be quantified as a k^~ 
dependent factor, C F {k 3 ) <x kf 2 . Since extragalactic 
point sources are expected to be the most significant as¬ 
trophysical foreground contaminate we include this con¬ 
tribution in our analysis as a nuisance term. 


2 . 2 . Reference Experiments 

The selection of reference experiments is critical to the 
outcome of the forecasting calculations since the uncer¬ 
tainties in redshifted 21 cm power spectrum measure¬ 
ments are coupled very closely to instrument design. The 
first generation of low-frequency radio arrays developed 
to perform redshifted 21 cm measurements includes three 
experiments: the Low Frequency Array 3 (LOFAR) being 
constructed in the Netherlands, the Primeval Structure 
Telescope 4 (PAST) under development in northwestern 
China, and the MWA-LFD 5 in western Australia. The 
designs of these instruments reflect different priorities 
and different approaches, but they all result in similar 
predicted sensitivities for measurements of the redshifted 
21 cm power spectrum. For example, although LOFAR 
will have the largest total collecting area (by a factor 
of ~ 10 ), its planned antenna distribution scheme is bet¬ 
ter suited for imaging applications than power spectrum 
measurements and its sensitivity is consequently less¬ 
ened. The design of the MWA-LFD, on the other hand, 
has been optimized for power spectrum measurements 
and exploits the wide fields-of-view that are inherent in 
dipole-based systems to overcome its more limited col¬ 
lecting area. 

Because of its focus on measurements of the redshifted 
21 cm power spectrum, we have chosen to base the design 
of the reference experiments in our calculations on the 
MWA-LFD. We will consider two experiments. The first 
is identical to the fiducial MWA-LFD design described 

3 http://www.lofar.org 

4 IPen et all f200l) 

5 http://haystack.mit.edu/ast/arrays/mwa 



Fig. 1.— Uncertainties (stepped lines) in spherically averaged 
power spectrum measurements due to thermal and cosmic sample 
variance for a single Az ~ 0.5 (8 MHz) region at redshift z = S 
observed by the MWA and MWA5000 reference experiments rel¬ 
ative to the fiducial, redshifted 21 cm power spectrum (dashed 
line). The thermal uncertain (thin gray lines) is overshadowed by 
uncertainty due to cosmic sample variance at large scales (small 
k). The combined uncertainty is shown as black lines. Due to the 
limited redshift range over which cosmic evolution can be neglected 
(Az ~ 0.5), the measurements in the light-gray region to the left 
of the vertical line are constrained only by angular fluctuations 
and will thus be most affected by astrophysical foreground con¬ 
tamination. The predicted constraints on the cosmological model 
parameters in Figure [2] and Tables [5] and [6] were derived using the 
full three-dimensional power spectrum, and thus, this plot of the 
uncertainty in the spherically averaged measurement is best used 
as a guide to illustrate which scale sizes will be most important in 
the constraints. 

in iBowman et all (|2006t Their Section 2). The array 
design consists of N = 500 antennas distributed within a 
D = 1500 m diameter circle. The density of antennas as 
a function of radius is taken to go as ~ r , but capped 
at a maximum density of one antenna per 18 m 2 . The 
antenna response is approximated by 

W{6) = cos 2 (§0/0), 9 < 0 (3) 

where 0 is proportional to wavelength and spans from 
31 to 38° over the range of interest. The angular resolu¬ 
tion of the array is given by A /D and the total collecting 
area by N dA , where dA is the collecting area of each 
antenna and scales like dA = 16(A 2 /4) for A < 2.1 m 
and is capped for longer wavelengths. Finally, the full 
bandwidth of the instrument is B = 32 MHz and the fre¬ 
quency resolution is 8 kHz. All of the fiducial properties 
are summarized in Tables [T] and [2] 

For the analysis in this paper, we define the observa¬ 
tion to be of a single field with 1000 hours integration 
during the most favorable circumstances. Additionally, 
we set the frequency coverage to 125 < / < 157 MHz, 
which spans 10 > z > 8 , and treat the bandwidth as 
four consecutive 8 MHz regions, each of which spans ap¬ 
proximately Az = 0.5. In principle, the observed 21 cm 
power spectrum is varying continuously with redshift due 
to cosmology-dependent effects. By dividing our mea¬ 
sured volume of space into thin regions in redshift we 
may ignore cosmic evolution within each region. 
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The second reference experiment is based on an ex¬ 
panded MWA -LF D configu ratio n similar to the one con¬ 
sidered by IWvithe et al.l (l2005lh although more con¬ 
densed. Dubbed the MWA5000, the expanded array 
contains 5000 antennas distributed over a 3000 m di¬ 
ameter region. The density of antennas remains ~ r -2 
and capped at one per 18 m 2 . The antenna response, 
instrument bandwidth, and frequency resolution are also 
unchanged (see Tables [l] and [2j • In both experiments 
all antenna elements are correlated to preserve the large 
field of view. 

Figure [T| illustrates the relative sensitivities of the 
two reference experiments by plotting their uncertain¬ 
ties in spherically averaged bins due to thermal noise 
and cosmic sample variance along with a fiducial, red- 
shifted 21 cm power spectrum. The two reference ex¬ 
periments share a common level of uncertainty at large 
spatial scales (small k ) due to cosmic sample variance 
since their fields of view are identical. At smaller spa¬ 
tial scales (k > 0 . 1 ), however, the larger collecting area 
of the MWA5000 reduces its thermal uncertainty in a 
given integration time compared to that of the MWA. 
The MWA5000 is nearly an order of magnitude more 
sensitivity than the MWA at these scales. 

3. METHOD 

The Fisher information matrix provides a convenient 
method for translating uncertainties in power spectrum 
measurements to cons traints on cosmological parameters 
(|Tegmark et al.l[l997f) . The minimum errors are calcu¬ 
lated using the Fisher matrix: 

p y- 1 dP(k t ) dP(k t ) 

ab 4" CT ( k *) 2 d Pa d Pb ’ U 

where er 2 = <j 2 n + cry is the combined uncertainty per 
independent cell due to thermal noise and cosmic sample 
variance, P is the measured power in a cell given by 
Equation^ and p is the set of model parameters. Taking 
the square roots of the diagonal elements of the inverse 
of Fisher matrix gives the errors. 

The model we use for the observed power spectrum 
represents the measurement after substantial data re¬ 
duction has occurred, including instrumental calibration 
and astrophysical foreground subtraction. At this stage, 
the observed power spectrum is parameterized based on 
three contributions: 

P(k) = P HI (k) + P F (k 3 ) + Pn, (5) 

where the first term is the redshifted 21 cm contribution, 
the second is the primary residual astrophysical (extra- 
galactic point source) foreground contribution discussed 
in Section 2.1, and the third is the residual white thermal 
noise contribution (also discussed in Section 2.1). 

For our assumption that reionization has heated the 
neutral hydrogen in the IGM but not yet produced sig¬ 
nificant features in the emission patterns of the gas by 
the target redshifts, the redshifted 21 cm contribution 
to the model power spectrum follows the matter power 
spectrum and can be expressed at a given redshift as 

PHiik) = Cj y Pss{k |A, !1m, fib, h,n s ,a s ), ( 6 ) 

where Cj y is the strength of the 21 cm emission from 
mean-density neutral hydrogen gas in the IGM and 


Pbb(k) is the matter power spectrum. At the target red- 
shifts, fluctuations in the matter density are expected 
generally to be in the linear regime over the range of 
spatial scales sampled. For redshift 2 = 8 , the pertur¬ 
bations at the smallest scales constrained by the mea¬ 
surements (k ~ 1 Mpc -1 ) will be approximately 40%; 
and at the largest scales (k ~ 0.01 Mpc -1 ), the pertur¬ 
bations will be less than 1%. In contrast to measure¬ 
ments of the matter power spectrum through large scale 
structure surveys of galaxies, redshifted 21 cm measure¬ 
ments probe directly the baryon density perturbations 
in the IGM, which follow closely the dark matter density 
perturbations. The effect of small deviations from the 
purely linear regime can be modeled easily in this case 
and we omit, therefore, this complication in our analysis. 
We use six parameters, p = {A, CIm, fib, h, n s , a s }, to de¬ 
scribe the matter power spectrum. These parameters are 
summarized in Table [3] along with their fiducial values. 
Additionally, we constrain Qk = 0 and = 1 — Hm- 
To compu te t h e matter pow er spe ctrum, we use CMB- 
FAST (jSeliak fe Zaldarriagalll996l ) and do not include 
velocity distortions or deviations from linear dynamics 
in the model. 

As described in Section 2.1, the primary residual fore¬ 
ground contribution is a function of k^ and is expressed 
as 

P F (k 3 ) = Fx c norm (, ( 7 ) 

where F is an amplitude scale factor of order unity and 
Cnorm, is a normalization constant. The thermal offset 
contribution does not vary with k and is defined as 

Pn — N x c norm , ( 8 ) 

where N is also an amplitude scale factor of order unity. 
For the Fisher matrix analysis, both P F and Pn are nor¬ 
malized by setting k n0 rm = 0.01 Mpc -1 and evaluating 
Cnorm = PHi{k nor m ) at redshift 2 = 8 . Thus, the mag¬ 
nitudes of F and N give the amplitudes of the residual 
contributions approximately relative to the peak value 
of the redshifted 21 cm power spectrum. Adding F and 
N to the set of six cosmological parameters yields the 
eight parameters, p = {A, SIm, fib, h , n s , a s , F, N}, which 
constitute the complete set of free variables of our base 
model. 

4. RESULTS 

The results of performing the Fisher matrix calcula¬ 
tions for the two reference experiments described in Sec¬ 
tion [272] are summarized in Tabled which lists the fore¬ 
casted 1-cr uncertainties on the model parameters. These 
values indicate that observations with the MWA-LFD 
would not constrain the cosmological parameters signifi¬ 
cantly. The only exception to this finding is for a s , which 
could be constrained to 0 ± 0.04 under the favorable as¬ 
sumptions used in the model. The second reference ex¬ 
periment, MWA5000, does constrain reasonably all the 
model parameters. Table U shows that this reference 
experiment provides constraints at levels approximately 
eq uivalent to those from the first-year WMAP results 
(jS per gel et alJ120031 Their Table 10). 

The elements of the covariance matrices for each refer¬ 
ence experiment are given in Tables [5] and [Gl and addi¬ 
tional information about the forecasted constraints for 
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Fig. 2.— Marginalized elliptical error regions for pairs of model parameters for the MWA5000 reference experiment. The contours are for 
68% (solid) and 95% (dashed) likelihood. The first two columns are set apart to emphasize the distinction between the nuisance components 
in the model due to astrophysical foregrounds and the cosmological terms. The bounds of the plotted regions for the cosmological parameters 
are set at twice the uncertainties (zb 2-cr) reported in the WMAP first-year results for the respective parameters. Thus, if the dotted ellipse 
is visible in a particular plot, the constraints on the parameters would be an improvement over WMAP. It is evident that the MWA5000 
would do a relatively good job of constraining 9 m, n s , and a s , and a comparatively poor job of constraining and h. The inability of 
the MWA5000 to constrain the scalar amplitude, A, of the power spectrum is due to degeneracies with other parameters, in particular 
and h, but also n s . 


the MWA5000 is provided in Figure 0 which shows 
marginalized error ellipses for two-parameter combina¬ 
tions of the model parameters. The properties of these 
results are discussed in detail in the remainder of this 
section. 

4.1. Parameter Dependencies 

To facilitate a conceptual understanding of the origin 
of the constraints on the model parameters, it is instruc¬ 
tive to individually vary the cosmological parameters and 
compare the resulting redshifted 21 cm signals. On all 
but the largest length scales, the experimentally observed 
power spectrum is well approximated by 


where W is the integrated value of the observational 
window function and represents the signal strength for 
a given instrument, and d^k/dudvdr/ is the Jacobian for 
converting the units of the matter power spectrum to the 
obse rved units. The relation ships between u,v,r ] and k 
are dMorales fe Hewittll2004D : 

27T 

[ki,k 2 \ = D [u,v\ ( 10 ) 

2ttv 21 H 0 E(z) 

* 3 “ c(l + z)2 <U) 

where 

e(z) = y/n M { i + zf + n k (i + z) 2 + n A (i + s) 3 ^™), 

( 12 ) 


















































6 


x 10 8 



Fig. 3.— Redshifted 21 cm power spectrum for the fiducial 
cosmological model parameters, shown with four example varia¬ 
tions of the parameter values. The top and bottom panels show 
the same elements plotted in different units. The thick solid line 
is for the fiducial model with standard cosmological parameters 
(Q b = 0.05, n M = 0.30, h = 0.70, n s = 1.00, a s = 0). The 
other curves are produced by adjusting one cosmological param¬ 
eter from the fiducial: Q b + 10% (dotted), h + 10% (dash), n s 
+ 10% (dash dot), and a s =0.1 (thin solid). The dark shaded 
region indicates the 1-cr uncertainty due to cosmic variance within 
the observed field and the light shaded region (at large k) is the 
thermal uncertainty for the MWA5000 reference experiment. As 
in Figure [I] the measurements to the left of the vertical line are 
constrained by angular fluctuations only and will be most sensitive 
to astrophysical foreground contamination. 


and Dm{z) and z^i = 1420.4 MHz are the transverse 
comoving distance and the rest-frame frequency of the 
hyperfine line, respectively. 

In Equation [9] there are three separate contributions 
which may be affected by changing the cosmological pa¬ 
rameters: the matter power spectrum itself, the neutral 
hydrogen emission factor, and the coordinate mapping 
between k and it, v, 77. The mapping between k and 
u,v,7y produces two effects when measuring the power 
spectrum. First, the region in it, v and 17 occupied by 
the observed power spectrum depends on the cosmo¬ 
logical model. This effect is analogous to the observed 
CMB power spectrum shifting in t for different cosmo¬ 
logical models, but since the redshifted 21 signal is three- 
dimensional, there may be a difference in scaling be¬ 
tween the sky-plane and line-of-sight dimensions. This 
additional degree of freedo m gives ris e to the AP test 
( Alcock fe Paczvnskil Il979f) (see also ISeo fc Eisensteinl 
( 2003f) for an applicable discussion of constraining dark 
energy using baryon acoustic oscillations in power spec¬ 
tra from three-dimensional galaxy redshift surveys). 

The second effect of coordinate mapping in the ob¬ 
served power spectrum is the presence of a significant 
cosmology-dependent amplitude factor in the measured 
power spectrum due to the Jacobian contribution in 
Equation [9] Differentiating and combining Equations [lO] 


and [Til yields the Jacobian: 

d 3 k = (27r) 3 V2 iHqE{z) 
dudvdrj C (1 + zf D 2 M (z) 

The effect of the Jacobian amplitude factor is further 
compounded by the hydrogen emission factor, Cj y . This 
term is unique to redshifted 21 cm experiments and pro¬ 
vides a significant boost to the cosmology-dependent am¬ 
plitude factor. Assuming the spin temperature of the 
neutral hydrogen in the IGM is significantly warmer 
than the CMB, the hyd rogen emission factor goes as 
(|Morales fc Hewittl 12004 Their Equation 13 and Their 
Appendix A) 

Ufo/i _ . . 

Cjy a ~E(z ) Xm ’ ( 14 ) 


where Shi is the mean neutral fraction (taken to be unity 
in this paper). Using Equations [13] and [14] and recalling 
that Dm{z) oc /i , the amplitude of the observed power 
spectrum is now seen to be proportional to 


Phi oc 


nlh 5 

(1 + z) 2 E(z) 


Pss- 


(15) 


Equation [15] indicates that the amplitude of the red¬ 
shifted 21 cm fluctuations does contain useful informa¬ 
tion about the cosmological model (with the caveat that 
the HI spin temperature must be known). This is in 
contrast to measurements of the matter power spectrum 
by galaxy clustering surveys, where the connection be¬ 
tween the observed amplitude of the power spectrum 
and the underlying physical amplitude remains weakly 
understood. 

Figure [3] shows the spherically-binned redshifted 21 cm 
power spectrum signal for the fiducial model (thick line, 
n b = 0.05, n M = 0.30, h = 0.70, n s = 1.00, a s = 0) 
along with the expected signals obtained by varying in¬ 
dividual parameters. We can see that varying or h 
primarily results primarily in a simple scaling of the ob¬ 
served power spectrum, as expected from Equation fl5l 
Thus these parameters are highly degenerate with each 
other and with the amplitude factor A, as shown in Fig¬ 
ure [2] and Table [ 6 ] 

The primordial power spectrum slope n s and running 
of the spectral index a s do not enter into the overall 
amplitude factor in Equation [15] for the observed power 
spectrum, but instead affect the shape of the underly¬ 
ing matter power spectrum. Increasing n s lowers the 
power at small k and increases the the power at large k, 
whereas increasing a s boosts the power spectrum at both 
extremes. Because these parameters do not strongly af¬ 
fect the amplitude normalization, they are relatively in¬ 
dependent of A,n 6 , and h as shown in Figure [2] and Table 
[G] They are also independent of our foreground parame¬ 
ter, F, because the shape they introduce is different than 
the k 3 2 shape of the residual foreground contamination. 
Thus n s and a s are the best constrained cosmological pa¬ 
rameters for the planned redshifted 21 cm observations. 

The last class of observational effects due to the cos¬ 
mological model is from the coordinate mapping. Most 
significantly, for an increase in h , a radially inward shift 
of scales (see Equations fTOl and fill with Dm{z) oc /t _1 ) 
counteracts the amplitude increase from the Cj y and 
Jacobian factors over much of the observed range in k- 
space. Increasing 9 m also shifts scales inward (through 
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E(z ) and Dm{z)), but does so differently for the line-of- 
sight direction and the transverse directions. This AP ef¬ 
fect cannot be represented easily in Figured however, we 
have included it in all the constraints. The significance 
of the effect is limited since we take Cl a = 1 — CIm- A 10% 
increase in CIm yields an approximately 4% reduction in 
scale, but only about 1% distortion between the direc¬ 
tions. In general, coordinate based effects may be better 
constrained with observations targeting the quasar bub¬ 
bles during reionization since sharper spatial features in 
the power spectrum may emerge during reionization due 
to characteristic sizes of reionized regions and Stromgren 
spheres around quasars. 

4.2. Degeneracies in Amplitude 

The uncertainty in A is considerable. Table |4] shows 
that the MWA5000 only constrains A to within 45%. The 
origin of this uncertainty is explained by the amplitude 
factor discussed above and is clearly evident in Table El 
which displays the large degeneracies (as normalized co- 
variance factors with magnitudes close to unity) between 
A, Clb, h, and more unexpectedly, n s . The degeneracy 
with n s is due to the range of scales to which the exper¬ 
iments are sensitive, which is primarily after the pivot 
point for the primordial power spectrum spectral index 
(k = 0.05 Mpc -1 ). Therefore, changes in n s tend to ap¬ 
pear less similar to tilting to the power spectrum and 
more like changing the amplitude, although n s remains 
well constrained. Additionally, the covariance between A 
and Clb confirms that the baryon acoustic oscillations in 
the power spectrum are not contributing much informa¬ 
tion about the baryon density relative to the amplitude 
changes from the Cj y and Jacobian factors. 

4.3. Residual Foreground Contamination 

Mitigating astrophysical foreground contamination is 
an important issue for redshifted 21 cm experiments, 
whether attempting to constrain cosmological models or 
determine reionization processes and scenarios. To this 
end, the relationships between Phi , Pf, and Pm, as de¬ 
scribed by the covariance between our model parameters 
A, F, and N are very useful to study in some detail. 

The constraints on the amplitudes of the residual as¬ 
trophysical foreground and residual thermal offset terms 
are especially encouraging. Neither contribution is signif¬ 
icantly degenerate with changes to the redshifted 21 cm 
power spectrum. In particular, the long, diminishing tail 
of the redshifted 21 cm power spectrum provides many 
statistical samples which contribute little to knowledge of 
the cosmological parameters, yet provide good references 
to prevent degeneracy with changes in the magnitude of 
the residual thermal offset. Thus, Pm is constrained to 
better than 2% of the peak value of Phi for the MWA- 
LFD and to better than 0.01% for the MWA5000. 

As discussed in Section m the residual foreground 
contamination is expected to be a power-law along the 
line-of-sight direction and to have little structure in the 
transverse directions. This is distinct from the redshifted 
21 cm power spectrum with its generally spherical sym¬ 
metry, although geometrical effects due to the shape of 
the observed volume of space work to negate this distinc¬ 
tion. The observed region has dimensions of ~9000 by 
~9000 Mpc in the sky-plane, but only ~100 Mpc along 



Fig. 4. — Unmarginalized error ellipses for the amplitudes of the 
redshifted 21 cm power spectrum and primary residual foreground 
contamination, A and F, forecasted for the MWA-LFD. The thick 
curves are for a Fisher matrix calculation utilizing the full three- 
dimensional measured power spectrum, while the thin curves are 
for the same calculation using spherical bins in /c-space. The con¬ 
tours are for 68% and 95% likelihood. The behavior is similar for 
the MWA5000, although the errors are smaller. 

the line-of-sight dBowman et al.ll200f)[ ). Since all three di¬ 
mensions are sampled with an approximately equal num¬ 
ber of divisions, this flattened shape is inverted when 
transforming into fc-space, leading to a highly elongated 
data cube in the ^ 3 -direction. For the MWA-LFD, the 
maximum k in the sky-plane (^ 1 ,^ 2 ) is approximately 
0.05 Mpc -1 , while along k$ it is over 10 Mpc -1 . There¬ 
fore, much of the information in the measurement about 
the redshifted 21 cm power spectrum is actually com¬ 
ing from the line-of-sight direction. Furthermore, the 
lack of depth in the observed volume prevents the ex¬ 
periments from probing large scales in the line-of-sight 
direction and, thus, the peak of the power spectrum (at 
k ss 0.01) is not well constrained. This means that not 
only is most of the information coming from the line- 
of-sight direction, it is coming from the power law-like 
tail (~ k~ 3 ) of the power spectrum. These effects com¬ 
bine to make the redshifted 21 cm contribution and resid¬ 
ual foreground contamination less distinguishable in the 
measured power spectrum. 

Fortunately, the analysis indicates that the residual 
foreground component is well separated from the red¬ 
shifted 21 cm power spectrum and does not affect sub¬ 
stantially the ability to constrain the cosmological pa¬ 
rameters. From Table El and Figure El we see that the 
largest degeneracy is between F and CIm- Removing Pp 
from the model power spectrum, and thus F from the 
parameter set such that p ={A, CIm, Clb, h, n s , a Sl N}, 
would reduce the forecasted uncertainty on CIm by less 
than 30% in an observation with with MWA5000. The 
contributions from the residual thermal offset, Pm, are 
even less. 

In our analysis so far we have treated the full three- 
dimensional measured power spectrum. This has maxi- 





mized the effects of the symmetry differences between the 
21 cm power spectrum and primary residual foreground 
power spectrum. We might have instead considered a 
reduced one-dimensional power spectrum produced by 
averaging over spherical bins in fc-space. This approach 
is used commonly to reduce data in which the expected 
signal has approximately spherical symmetry. 

In Figure [4) we compare the full three-dimensional 
treatment with the spherically-binned approach by plot¬ 
ting unmarginalized error ellipses in the A-F plane of 
parameter space for both cases. It is clear that the 
binned method produces a significantly larger degener¬ 
acy between the magnitude of the residual foreground 
contribution and the amplitude of the 21 cm neutral hy¬ 
drogen power spectrum than the full three-dimensional 
treatment. 

4.4. Dark Energy Equation of State 

Determining the nature of dark energy has become 
an important goal in astrophysics and the applicabil¬ 
ity of future large-scale structure surveys to this topic 
has b een considered recently with generally favorable re¬ 
sults (ISeo fe Eisensteinll2003tlLindedl2003tlHu k. Haimanl 
1 200. 'ft IWang et al.l 2004 lAbdalla fc Rawlingl ~f2005l . 
With their similarities to large-scale structure surveys, 
would measurements of the redshifted 21 cm neutral hy¬ 
drogen power spectrum also improve knowledge of dark 
energy? 

Including the dark energy equation of state parame¬ 
ter, w, as a free variable in the Fisher matrix analysis, 
so that p = {A,riM,D,b,h,n s ,a s ,w, F, N}, yields little 
information about the nature of dark energy. Neither 
the MWA-LFD nor the MWA5000 are able to provide 
meaningful constraints on w without priors from other 
experiments. Furthermore, the constraints on Dm, D;,, 
and h are significantly weakened (by up to an order of 
magnitude), although information is retained about the 
primordial power spectrum. 

Although the MWA5000 reference experiment would 
be able to detect the baryon acoustic oscillations (see 
Figure[3]), the experiment targets emission from high red- 
shifts prior to the epoch of dark energy domination. The 
value of the experiment would be relevant only as an 
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in intermediate measurement betwee n the C MB and low 
redshift large-scale structure surveys (lBarkanal2006l ). but 
could be important if dark energy had unusual behavior 
at high redshifts. 

5. CONCLUSION 

The first generation of redshifted 21 cm experiments 
have been shown previously to have the potential to char¬ 
acterize the processes and history of reionization at red- 
shifts 6 < z < 12. In this effort, we have considered a 
scenario in which reionization is not the dominant con¬ 
tribution to the measured power spectrum during this 
period. Our findings demonstrate that the initial exper¬ 
iments would not contribute substantially to knowledge 
of the underlying cosmology under these circumstances. 
While they may provide some use to cosmological stud¬ 
ies, the primary science return of the first generation of 
redshifted 21 cm experiments will be regarding the as¬ 
trophysics of reionization. 

Under the optimistic assumptions used in this pa¬ 
per, the second-generation experiments, similar to the 
MWA5000 reference experiment, will be better suited 
to constrain Dm and the primordial power spectrum 
through n s and a s . Even for the spectral index proper¬ 
ties it is difficult to envision a scenario where redshifted 
21 cm observations alone could provide an unambiguous 
constraint on cosmological information, however, since 
the signal could easily be confused if the universe had a 
small amount of reionization structure on small scales or 
spin temperature fluctuations in the IGM. 

We have also considered the covariance between the 
redshifted 21 cm power spectrum due to a fully neu¬ 
tral IGM and the primary residual foregrounds of low- 
frequency radio observations. The signal and contami¬ 
nation components were easily distinguished in the full 
three-dimensional power spectrum analysis. 
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TABLE 1 

Fiducial Observation Parameters 


Parameter 

MWA-LFD (MWA5000) 

Array configuration, p(r) (m —2 ) 

r>_> r~2 

Array diameter, D (m) 

1500 (3000) 

Bandwidth, B (MHz) 

32 

Frequency resolution (kHz) 

8 

Number of antennas, N 

500 (5000) 


Note. — Array parameters for the MWA-LFD and 
MWA5000 reference experiments. 
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TABLE 2 

Redshift Dependent Parameters 



z = 8 

z — 10 

Angular resolution (°) 

0.073 (0.036) 

0.089 (0.044) 

Antenna collecting area, dA (m 2 ) 

14 

18 

Antenna response scale, © (°) 

31 

38 

Frequency (MHz) 

158 

129 

System temperature, T svs (K) 

440 

690 


Note. — Characteristics of the fiducial observation that depend 
on frequency, and thus on redshift, are listed along with their values 
at the z = 8 and z = 10 edges of the observation redshift range. In 
the calculations, the values of the fiducial parameters at intermedi¬ 
ate redshifts were linearly interpolated from the end points (except 
for ©, see text). For the angular resolution, the MWA-LFD and 
MWA5000 have different properties due to the larger size of the 
MWA5000 (whose values are listed parenthetically). 



12 


TABLE 3 

Model Power Spectrum Parametrization 


Parameter 

Symbol 

Fiducial Value 

Normalization (at k = 0.05 Mpc -1 ) 

A 

1.00 

Matter density 

Qm 

0.30 

Baryon density 

f!;, 

0.05 

Hubble constant 

h 

0.70 

Scalar spectral index (at k = 0.05 Mpc -1 ) 

n s 

1.00 

Running index slope (at k = 0.05 Mpc -1 ) 

a s = dn s /d\nk 

0.00 

Residual foreground amplitude 

F 

0.00 

Residual thermal offset 

N 

0.00 


Note. — Eight parameters used to describe the measured power spectrum, 
and their values for the fiducial model. The first six rows give the basic param¬ 
eters used to describe the matter power spectrum and the last two rows give 
the parameters to quantify the residual astrophysical foreground and residual 
thermal offset terms. 
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TABLE 4 

Forecasted Uncertainties for the Model Parameters 



A 

F 

D 

Om 

Slj, 

h 

n s 




1.00 

0.00 

0.00 

0.30 

0.05 

0.70 

1.00 

0.00 

MWA-LFD 

6.69 

0.09 

2e-3 

0.21 

0.07 

0.88 

0.2 

0.04 

MWA5000 

0.45 

2e-3 

le-4 

0.01 

6e-3 

0.05 

0.02 

4e-3 

WMAP 

0.09 



0.04 

4e-3 

0.04 

0.03 

0.02 

SDSS 

0.10 

— 

— 

0.05 

2e-3 

0.04 

0.03 

— 


Note. — Forecasted 1 -c uncertainties for the model parameters. 
The first row gives the fiducial values for reference. The MWA-LFD 
reference experiment does not significantly constrain cosmological 
parameters, while the MWA5000 does a reasonably good of job, 
particularly for and the primordial power spectrum spectral 
index descriptors, n s and a s . For comparison, the last two rows 
display approximate constraints for similar model s from WMAP 
dSpergel et al.l 120031 , Their Table 10) and SDSS dTegmark et al.l 
l2004 Their Table 3). 


TABLE 5 

Covariance Matrix for MWA-LFD 



A 

F 

D 

FIm 

fib 

h 

n s cxs 

A 


1.00 







F 

-0.09 

1.00 






D 

-0.09 

0.01 

1.00 





flM 

-0.04 

-0.72 

-0.04 

1.00 




ft b 

-0.87 

-0.24 

0.04 

0.49 

1.00 



h 

-0.97 

0.24 

0.11 

-0.20 

0.72 

1.00 


n s 

0.95 

0.04 

-0.15 

-0.12 

-0.82 

-0.92 

1.00 


-0.41 

-0.37 

-0.15 

0.63 

0.69 

0.23 

-0.32 1.00 


Note. — Elements of the covariance matrix for the MWA- 
LFD calculated by taking the inverse of the Fisher matrix, F. 
Each element has been normalized according to c a b/y/c a aCbbi 
where the diagonal elements are the squares of the 1-cr uncer¬ 
tainties given in Table [4] 


TABLE 6 

Covariance Matrix for MWA5000 



A 

F 

D 

flM 

fib 

h 

n s cxs 

A 


1.00 







F 

0.20 

1.00 






D 

-0.07 

0.00 

1.00 





flM 

-0.40 

-0.60 

-0.02 

1.00 




fl b 

-0.93 

-0.35 

0.05 

0.66 

1.00 



h 

-0.97 

-0.08 

0.08 

0.19 

0.81 

1.00 


n s 

0.93 

0.27 

-0.07 

-0.46 

-0.83 

-0.92 

1.00 

ol s 

-0.48 

-0.28 

-0.20 

0.53 

0.50 

0.42 

-0.62 1.00 


Note. — Same as Table 3, but for the MWA5000. 






























